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Abstract 


A numerical simulation of .diode, using ^finite diflference method , is 
carried out for the turn ON and turn OFF transient conditions. The method is 
applicable to long-base and short-base diode for any type of doping profile. For 
a specified input waveform of current or voltage as a function of time, the pro- 
gram yields terminal currents and various quantities of interest in the interior 
of the device such as minority carrier density, different current components, 
electric field, electrostatic potential etc, as a function of time and position. 
The effect of changing external voltage, load resistance and the lifetime of the 
carriers on the turn OFF transient is investigated and comparison with ana- 
lytical results is made. The comparison between simulated and SPICE results 
for the turn OFF transient is also presented. 
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Chapter 1 


Introduction 


Most of the semiconductor devices contain at least one junction between 
p-type and n-type semiconductor materials, pn junctions are fundamental 
to the performance of many electronic circuits. The pn junction diode is a 
two terminal device which can perform various functions: (i) it is capable of 
rectifying low' frequency alternating signals, (ii) because of non ohmic nature 
of the I-V curve a pn diode can be used as a variable resistor (varistor), (iii) 
since the depletion layer capacitance of such pn junction can be made to vary 
as some suitable function of reverse bias applied to them, the diodes have been 
used as a variable capacitors (varactors), (iv) due to non linear nature of I-V 
characteristic, pn junction diodes are used in demodulation of communication 
signals, (v) if the terminals of such junction are left open and the junction is 
irradiated externally, carriers are generated, and according to their polarity 
they are collected by appropriate ends of the depletion layer, this causes an 
open circuit voltage to develop across the pn junction. Hence they can be used 
as solar cells, (vi) direct bandgap semiconductors can be used as LEDs and 
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semiconductor lasers. 


An important feature of a semiconductor diode is the similarity of its I-V 
characteristics with an ideal switch, hence, diodes are very important element 
in pulse electronics (mainly as a switching element in computers). High speed 
pulse diodes are also very important in radar and other branches of high and 
super-high frequency electronics. The rapid advancement in computer tech- 
nology led to the development of new construction principle of logic, memory 
and other computer elements, which call for switching of large current. It 
therefore becomes necessary to create high speed diodes capable of switching 
on currents in a very short period of time. Hence, pulsed diodes with high 
frequency and power are required. 

The development of new semiconductor technologies and novel semicon- 
ductor device structures has traditionally been guided by experimental ap- 
proaches. Starting from an established process sequence, fabrication steps are 
changed together with geometrical dimensions (feature size). The modified 
process is then realized by fabricating several lots. The finished devices are 
tested to ensure that their performance conforms to the initial design speci- 
fications. This approach usually includes several iterations of the process to 
testing loop. With the advent of increasingly complex integrated circuits, the 
traditional empirical approach becomes expensive and time consuming. 

In order to characterize the devices in a reasonable way, numerical methods 
emerge out to be an excellent solution. In the numerical methods, the actual 
device behaviour is simulated, taking into account semiconductor physics in 
detail. This is, of course, an intensive procedure in terms of the CPU time. 
However, it yields very useful information. 
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Numerical simulation of a semiconductor device is carried out by solving 
three non-linear partial differential equations, i.e. Poission’s equation and 
two continuity equations, one for electrons and the other for holes, inside the 
geometry of the device. Modelling based on partial differential equations, 
w'hich describe different regions of semiconductor devices in a unified manner 
was suggested by Gummel[l] for the one dimensional bipolar transistor. This 
approach was further developed and applied to pn junctions by De Mari[2], 
and to IMPATT diode by Scharfetter and Gummel [3]. 

The present thesis work deals with the simulation of a one dimensional 
silicon pn junction diode for the switching conditions. In Chapter 2, basic 
semiconductor equations, their discretization and the boundary conditions for 
the different types of excitations are described. The results obtained from the 
simulation and the analytical expressions are described in Chapter 3. In the 
last chapter (Chapter 4) conclusion and scope for future work are given. 



Chapter 2 


Basic semiconductor Equations 
and Device Simulation 


2.1 Basic semiconductor equations 


To determine the behavior of an arbitrary semiconductor structure under vari- 
ous operating conditions, a mathematical model is required. The semiconduc- 
tor devices operate under non-equilibrium conditions in which electron and 
hole current flow through the device. The carrier current at any point in the 
device can be calculated by the electric field and the carrier concentration gra- 
dient at that point. Hence, to calculate the quantities of interest we require 
appropriate mathematical expressions. 

The basic semiconductor equations are used to describe the static and dy- 
namic beha\ior of carriers in the semiconductor under the influence of external 
fields which causes the deviation from the thermal equilibrium conditions. The 
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three basic semiconductor equations [4] are given as; 


1. Poisson’s equation 


(2.1) 

e 

where, ip = electrostatic potential, p = space charge density = g{p — n + — 

^^a)i ~ ionized donor density, N~ = ionized acceptor density. 

The Poisson’s equation relates the potential to the charge density and thus 
enables us to determine the electric field for a given charge distribution. 


2. Continuity equation for electrons 

V.(I?nVn - pnuVip) - R{n,p) = dn/dt (2.2) 

3. Continuity equation for holes 


V.(jDpVj9 + PpP^ip) - R(ji,p) = dp/dt 


(2.3) 


where Dn = diffusion constant for electron, Dp = diffusion constant for holes. 


R = recombination rate which is given by the Shockley-Read-Hall formula 


R 


np - 


(2.4) 


Tno{p + Pt) + rpo{n + nt) 

where Ut = electron trap density, pt = hole trap density, Xno = electron lifetime 
Tpo = hole lifetime. 


The continuity equation expresses the principle of particle conservation by 
relating the net rate of change of charge carriers in a volume to the process 
that introduce and remove the carriers. 


The Poisson’s Equation is elliptic and the continuity equation is parabolic 
in nature. To analyze the behaviour of the device, above system of partial 
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differential equations is to be solved, however, in general, these equations can- 
not be solved analytically. Hence numerical methods are used to solve these 
equations. There are three tasks involved to solve these equations which are 
as follows: 

1. The domain, i.e., the simulation geometry of the device is divided into 
finite number of partitions (grid points), -where the solution can be approxi- 
mated with the desired accuracy. 

2. The differential equations are approximated in each of the domain by 
algebraic equations w'hich involve only values of the continuous dependent 
variables at discrete points in the domain, and the knowledge of the structure 
of the chosen function -w'hich approximate the dependent variables within each 
of the sub domain. Hence, a large number of nonlinear algebraic equations are 
obtained with unknowns comprised of approximations of continuous dependent 
variables at discrete points. 

3. The system of nonlinear equations is solved using the Newton Raphson 
method. 


2.2 Scaling 


The first step before the analysis of the semiconductor equations is proper 
scaling because the dependent variables (^,n,p) in the basic semiconductor 
equations have large difference in magnitude. After scaling all the variables 
involved in the semiconductor equations have nearly same range of variation. 
In the present work the scaling factors used are given in Table 2.1. 
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quantity 

symbol 

value 

X 

Xo 

max(| X 1) 

tp 

V’o 

kT/q 

n,p,C 

Co 

max| C{x) 1 

Dn,Dp 

Do 

max(D„,Dp) 

R 

Ro 

b 

o 

OlO 

t 

h 

Xq/Dq 


Table 2.1; Scaling factors 


After scaling the basic semiconductor equations becomes; 


\\V^-{h-p-C) = 0 

(2.5) 

- - dh 

V.(r>„Vn - fMnnVi>) - R{n,p) = 

(2.6) 

~ ~ dp 

V.(DpVp + PppVi)) - R{n,p) = 

(2.7) 


2.3 Discretization of basic semiconductor equa- 
tions 


There are several methods of deriving algebraic equations from the differential 
equations. In the present work, we are using the finite difference method. The 
whole analysis is carried out in one dimension. 

In the finite difference method, the domain in which the solution is to be 
obtained is partioned into subregions. In one dimensional case, the boundary 
of the domain is two points A and B, as shown in the Fig 2.1 . Between A and 
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Figure 2.1: Discretization scheme 


B is the interior of the device, where the solution is required. At each node 
the values of the unknowns are represented b)' t/’i, 'n-h Pi- For each node we 
have three unknowns and three equations and also at the boundary, we have 
appropriate boundary conditions for each unknown. Hence for N nodes we will 
get 3N equations. 

First we will derive the difference equations for the static semiconductor 
equations, which are obtained if the boundary conditions for the electrostatic 
potential are time invariant. The partial derivatives of electron and hole con- 
centration with respect to time vanish in this case. Hence the basic semicon- 
ductor equations becomes: 


A^V 2 ^-(n-p-C ')=0 

( 2 . 8 ) 

V.(D„Vn - - R{n,p) = 0 

(2.9) 

V.iDpVp + pppViP) - R{n,p)=0 

( 2 . 10 ) 


To derive the difference equations, the differential operators are directly 
replaced by the difference operators. For the Poisson’s equation, it is assumed 
that the electrostatic potential varies linearly between two adjacent points, as 
shown in Fig 2.2. Thus we can write 
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Figure 2.2: Finite difference scheme 


dx^ 


drp/dx |t+i /2 -dipidx |i-i /2 
2 


( 2 . 12 ) 


For continuity equations, the assumption of n and p varying linearly be- 
tween two grid points is not valid, because the percentage variation of n and 
p between the two grid points is much more than that of the electrostatic po- 
tential ijj. Hence, the discretization of the continuity equation is carried out 
using the Scharfetter and Gummel approximation [3]. In this approximation, 
it is assumed that the variation of n and p between two adjacent grid points 
is exponential in nature. The discretized continuity equation for electrons is 
given by [4] 

^ rB{i}i+i - ipi)ni+i jB(^i - ipi+i)ni _ 

hi hi hi-i 

I ~ '^ . ht ±)R{ri,p) = 0 

h,_i 2 


and the continuity equation for holes is given by [4] 

hi hi hi-i 
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_ ^ ^ytb zl)R[n,p) = 0 

(2.14) 

where B is Bernoulli function and is defined as [4] 

B(x) = ^ (2.15) 

2.4 Discretized equation for the transient con- 
dition 


In practical circuits, the external excitation varies with time, hence, elec- 
trostatic potential across the diode becomes time dependent and the partial 
derivatives of the carrier concentration with respect to time does not vanish. 
The discretization of these equation is obtained using fully backward time dif- 
ferencing (Backward Euler method) scheme. The discretized equations then 
becomes 

Fi (V>m+1 , ^+1 . Pm+l) = 0 (2-16) 

n„+i - = 0 (2.17) 

dm 

Pm+lZ lUL _ F 3 (V’m+l,nm+l,Pm+l) = 0 (2-18) 

dm 

where dm = tm+i -tm,tm and are the two consecutive time steps, 
nm,Pm and JV+iiPm+i are the concentrations of electrons and holes at time 
instant tm and Wi as shown in Fig 2.3, Fi is Poisson equation, Fa and Fg are 
the continuity equations for electrons and holes respectively. If we substitute 
the spatial operators in Eq 2.16-2.18, we get completely discretized equations 
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:entral library 
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Figure 2.3: Change in ip, n and p with time 


as 


\2r.(. hj^i + hi , ^hi-i + hi ^ 

^ \yi-l,m+l nr Wi,Tn+l\ ol oi, ) 


+^i+l,m+l 


/ij—i H" /ij 

2h, 


2/ii_i ' 2h, 

] “ Pi,m+1 ^i) ~ 0 


(2.19) 


ft' ft 

1pi-l,m+l) \i'^ + B{1p^,m+l - tpi+l,m+l)^" * + “ ) 


2hi-i 


2hi 




n. 


+«t+l,m+l-B(l/’t+l,Tn+l “ ~2/v ^ ~ ) 


( 2 . 20 ) 


^ I 

Dp\pi-i^rn+lB{lpi,m+l ~ 4’i-l,Tn+l) ^ ^ Pi, m+l (-6(^1-1,171+1 ' 

, .hi-i + hi , , hj-i + hi ^ 1 

^i,m+l) JTi! h 6(l/>i41,in+l ^i,m+l) nz. j ) 


2hi-i 


2K 


*'771 

+Pt+l,m+l6(^,,m+l - V'i+l.m+l) * 2/li ” (-^.m+l “ ^) 


( 2 . 21 ) 


2.5 Boundary conditions 

The above discretized equations cannot give the values of the electrostatic 
potential and carrier concentration at the boundary of the device. Hence, the 
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(2.27) 



(a) (b) (c) 

Figure 2.5: Different types of excitations. 


For a purely voltage controlled source, Fig 2.5a, if we apply and Vp on 
the left and right hand side boundary respectively, i.e, VappUed = yp -Vn, the 
electrostatic potential at the left hand boundary (V'l) will change to 

<pL = V„+,pl (2.28) 

and on right hand side boundary electrostatic potential {iPr) will become 

^R = Vp- (2.29) 

where and ipR are the values of the electrostatic potential at equilibrium, 

and are given by 

i,l = Vr!n(— ). (2-30) 

Hi 

and 

4 = W— ) (2-31) 

Tii 

For a purely current controlled contact, as shown in Fig 2.5b, the value of the 
electrostatic potential at the right hand boundarj' can be calculated using 

+ (2.32) 
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(2.27) 



(a) (b) (c) 

Figure 2.5: Different types of excitations. 


For a purely voltage controlled source, Fig 2.5a, if we apply Vn and Vp on 
the left and right hand side boundary respectively, i.e, VappUed = Vp — Vn, the 
electrostatic potential at the left hand boundary (V’l) will change to 

V’i = v; + 1/.S ( 2 . 28 ) 

and on right hand side boundar}' electrostatic potential (iPr) will become 

= (2.29) 

where and are the values of the electrostatic potential at equilibrium, 
and are given by 

= VHn(— ), (2.30) 

rii 

and 

= Vrlni^) (2.31) 

For a purely current controlled contact, as shown in Fig 2.5b, the value of the 
electrostatic potential at the right hand boundary can be calculated using 

j{Jn + Jp)-dA-I,{t) = 0 (2.32) 
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where A=cross sectional area of the device, is the current which is forced to 
flow through the device, Jn and Jp are the electron and hole current densities 
respectively, and are given by [4] 


Jn = D 


n 


- ‘ipR-i)nR 

hR~i 


(2.33) 


and 


Jp = Dp 


B{lpR-l - ^I}r)Pr - B{i}R - i>R-l)'pR-l 
hR-1 


(2.34) 


After putting the values of Jn and Jp in Eq 2.32, the electrostatic potential 
at the right hand boundary can be calculated. The potential at the left hand 
boundary (^l) is kept same the as the equilibrium value, which is given by 


ipL = VrlnipLlrii). 


(2.35) 


When the diode is used in a circuit having a load resistance as shown in 
Fig 2.5c, the boundary conditions for the electrostatic potential are obtained 
using KVL around the loop, given by 

V,-I R-Vd = 0 (2.36) 

where Vb is the voltage across the diode, and is given by 

VjD=1pR-1pL + Hi) (2.37) 

where Hi is the barrier potential and is given by 

H.=-0l + 4- (2.38) 

After the boundary conditions are established, the nonlinear equations are 
solved by the Newton-Raphson method [6]. 
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where A— cross sectional area of the dence, is the current w'hich is forced to 
flow through the device, and Jp are the electron and hole current densities 
respectively, and are given by [4] 


T _ n ~ — ipR-i)nR 



"•R-1 


(2.33) 


and 


A 


Bji^R-i - ^r)pr -BjjiR- 
h-R-i 


(2.34) 


After putting the values of and Jp in Eq 2.32, the electrostatic potential 
at the right hand boundary can be calculated. The potential at the left hand 
boundary {ipi) is kept same the as the equilibrium value, which is given by 


= VTlnipilni). (2.35) 


When the diode is used in a circuit having a load resistance as shown in 
Fig 2.5c, the boundary conditions for the electrostatic potential are obtained 
using KVL around the loop, given by 

V,-IR-Vp = 0 (2.36) 

where Vd is the voltage across the diode, and is given by 

VD = ^R-i>L + Vu, (2.37) 

where Vti is the barrier potential and is given by 

+ (2.38) 

After the boundary conditions are ^ablish^, the nonlinear equations are 

solved by the Newton-Raphson method [6]. 
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Chapter 3 


Results and discussion 


The simulation method as described in Chapter 2 can give the results for 
any arbitrary doping profile and external excitation. In the present work we 
have chosen abrupt doping profile because most of the analytical solutions are 
available for this type of profile. A sample device structure is shown in Fig 

3.1. 


- 

= 10 sec 
10 sec 


-3 


Nj 5 = 10 cm 

Na= 10 cm 

300Pm^ — ► 




Figure 3.1: Sample 


The length of the p side and n side of the diode are chosen ~ bin and ~ 5Lp 


15 




respectively so that the minority carriers should recombine before reaching the 
diode boundary. 


3.1 Turn on transient 


The phenomena occurring inside the device for the small current step and 
the high current step are quit different, which are explained in the following 
sections. The turn ON excitation (current or voltage) is showm in the Fig 3.2. 
It is defined as 


uit) = 


0 if t < 0 


F t>0 


F 


0 time 

Figure 3.2: Turn ON step. 


3.1.1 Small current step 


A small current step is recognized as the excitation causing low level in- 
jection on either side of the junction. We have chosen 


J(t) = I 


0 

10 Amp cm""^ 


if t < 0 

<> 0 
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Figure 3.3: Variation of diode voltage with time for small current step. 

For the above given current the injection level is of the order of 10^^ cm~^ 
(Fig 3.7, in all the subsequent figures ti = 1 ns, where ti is the time after 
the rising edge and each subsequent curv'e is At ns apart , where At = 1 ns), 
hence it is verified that the injection level is low. The voltage response for 
this excitation is shown in the Fig 3.3. The voltage waveform does not show 
an abrupt rise due to the charging of the transition capacitance [7]. The dif- 
fusion capacitance has insignificant influence on the transient response at low 
level injection. In the initial part, the transition region capacitance is charged 
by the displacement current because initially the displacement current is vary 
high in the transition region. In Fig 3.4 the displacement current is showm 
as a function of position and time. It has significant values near the junction 
because the variation in the electric field with time is significant only in this 
region. At later times, electron and hole currents become dominant compo- 
nents and the displacement current charging the transition region capacitance 
gradually decreases to zero. The variation of the electron and hole current in 
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the neutral p region are shown in Figs, 3.5 and Fig 3.6 respectively. 
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Figure 3.5: Electron current density vs distance for small current step. 
Prom Figs 3.4-3.6, it can be seen that initially the transition region is 
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Figure 3.6: Hole current density vs distance for small current step. 



300 310 320 330 

Distance (pm) 


Figure 3.7: Minority carrier (electrons on p side) concentration vs distance for 
small current step. 




Electric field ( Volts / cm ) 


800 








Figure 3.10; Electrostatic potential vs distance for small current step. 


dominated by the flow of displacement current. The electron and hole current 
components are very less in this region due to the little change of electric field 
and the mobile carrier densities from the initial equilibrium value (minority 
carriers on the p side are shown in Fig 3.7) while the variation of the electric 
field with time is large as shown in Fig 3.8. The variation of electric field with 
distance at different times is shown in the Fig 3.9. The decrease in the electric 
field in the transition region is related to the change of junction voltage. The 
electric field is maximum at the junction due to the maximum variation of 
electrostatic potential at that point. The electrostatic potential with respect 
to distance curve is shown in Fig 3.10. We see that due to some finite resistance 
of the neutral region the electrostatic potential curve is not constant in the 
neutral region. 
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3.1.2 High current step 


We have taken the following as a high current step 


J{t) = 


0 

500 Amp cm“^ 


if t < 0 
t > 0 


When it is applied to the device, initially in thermal equilibrium, the final 
state is high injection level. The minority carrier concentration is ~ 0.6. ( 
as shown in the Fig 3.16), which confirms high level injection. The variation of 
the diode voltage with time is shown in Fig 3.11. The initial part is the result 
of combination of two basic effects : the buildup of the ohmic voltage drop in 
the neutral regions and the charging of the transition region capacitance. In 
the high level injection the motion of the carriers is ruled by both drift and 
diffusion phenomena, as a result of conductivity modulation. The displacement 
currents at various instants of time in the transition region and the neutral 
region respectively are shown in Figs. 3.12 and 3.13. 

We see that this is different from the small current excitation. Figure 3.12 
shows the displacement current in the transition region. Due to the very small 
width of the transition region, the displacement current is nowhere constant 
as opposed to Fig 3.4 (where it is constant in the depletion region) for small 
current step, where the depletion width is larger than the present case. Figure 
3.13 shows the variation of displacement current in the neutral region at dif- 
ferent time instants. The reversal of the direction of the displacement current 
flow is a consequence of decrease of voltage drop in the quasi neutral region 
because conductivity modulation is occurring in this region. The modulation 
initiates at the edge of the transition region and evolves gradually towards the 
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Figure 3.11: Diode voltage vs time for high current step. 



Figure 3.12: Displacement current density vs distance for high current step in 
the transition region. 
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Distance (nm) 

Figure 3.13: Displacement current density vs distance for high current step in 
the neutral region 



Figure 3.14: Electron current density vs distance for high current step. 
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Figure 3.15: Hole current density vs distance for high current step. 



300 340 380 420 

Distance (pm) 


Figure 3.16: Minority carrier (electrons on p side) concentration vs distance 
for high current step. 
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Figure 3.17: Electrostatic potential vs distance for high current step, 
external contacts [7]. 

Figures 3.14 and 3.15 show the electron and hole current distributions on 
p side. We see that the slopes near the edge of the transition region is higher 
than those for small current step (Figs 3.5 and 3.6). 

The variation of the electrostatic potential is shown in Fig. 3.17. Due 
to the large voltage drop in the neutral region the slope of the electrostatic 
potential is high. 

3.1.3 Voltage step 

A voltage step (0 to 0.3 Volts) is applied to the diode without any load 
resistance to see the variation of the minority carrier concentration. For a 
voltage step, the variation of the diode current is shown in Fig 3.18. We 
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see that initially the current rises very sharply, this is due to the transition 
region capacitance (it acts as short circuit). Hence, initially the displacement 
current is also significant (Fig 3.19). As time progresses, the transition region 
capacitance starts to charge and the displacement current decreases, the total 
current also reaches the steady state value. The electrostatic potential is shown 
in Fig 3.20. Its variation occurs in the neutral region near the junction. The 
excess minority carrier concentration is shown in Fig 3.21. The excess electron 
concentration near the junction changes from zero to npo[exp(Vj?/Vr) - !]• 
Since the electrons cannot diffuse into neutral p region just after the turn ON 
step, the slope of the excess electron concentration at the edge of the transition 
region at the p side {{6n/6x) |x=xp) will be large just after the turn ON step, 
hence this is also a reason for the initial high current through the diode, and 
slowly as the minority carrier concentration increases the slope decreases and 
thus the current decreases [8]. 



Figure 3.18: Variation of diode current vs time for voltage step. 
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igure 3.19: 


Displacement current density vs distance for voltage step. 



Figure 3.20: Electrostatic potential vs distance for voltage step. 
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Figure 3.21: Minority carrier concentration (electrons on p side) vs distance 
for voltage step. 

3.2 Turn off transient 


The turn OFF step is shown in Fig 3.22a. For the turn off transient 
we have taken the circuit shown in the Fig 2.5c (Chapter 2), for analyzing 
the storage and fall time delay of the current through the diode, by varying 
different circuit parameters. The diode chosen is n^p for comparison with the 
analytical expressions. 

The applied input voltage and the diode responses are shown in the Fig 
3.22. Initially when the applied voltage is 10 Volts, the diode is forward biased 
and steady state current of 4.5 Amps is flowing (Fig 3.22c) through it. At time 
< 1 , the diode voltage is reversed to -5 Volts. The diode voltage reduces (Fig 
3.22b), but does not reverse. The current does not drop to zero, it reverses and 
remains at ~ V/Ri because the minority carriers cannot recover immediately. 
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At time #2 when the excess minority carriers in the immediate neighborhood 
of the junction have been swept back across the junction, the diode voltage 
starts to reverse and the magnitude of the diode current begins to decrease. 
Prom Fig 3.22d we see that from time ti to t 2 the rise in electric field is small 
because the transition width does not changes significantly, it starts rising only 
after time ta- 


The analytical formulae used for calculating storage delay time [9] and fall 
time [10] are as follows: 

1 


Ts = r„[erfc{ 


1 + ^r/ If 




erfi 


a/? 


(3.1) 


(3.2) 


where t„ is electron lifetime, Tg is storage delay time, rp is fall time, Ir-= 
reverse current and 1f= forward current. The effect of the variation of the 
load resistance on current is showm in Fig 3.23. We see that the storage 
time and the fall time delay remain constant. This is because the ratio of 
the forw'ard current (which is injecting the minority carriers) and the reverse 
current (which is withdrawing the minority carriers) remains constant. 


Figure 3.24 shows the variation of the current if the reverse bias is increased, 
keeping forward bias constant. We see that as the reverse bias is increased, 
both storage time and fall time decrease, which is due to the increase in the 
reverse current. As the reverse current increases, the minority carriers are 
swept back quickly. The variations of the storage time and the fall time delays 
are shown in the Figs. 3.25 and 3.26 respectively. When keeping reverse bias 
constant and forward bias changing, the current varies as shown in Fig.3.27. 
We see that the storage and the fall time delays are increasing. This is because 
of increased forward current and thereby increase in injected excess minority 
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carriers. The reverse current is constant, hence, more time is required to 
withdraw the excess minority carriers when the forward bias are larger. The 
increase of storage and fall time delays with increase in forward bias are shown 
in Figs 3.28 and 3.29 respectively. The analytical expression values of storage 
delay time agrees well with the simulated results but there is large difference 
betw'een the analytical results for fall time and the simulated results. This 
difference can be due to the fact that the change of the transition region width 
is not taken into account properly when deriving the analytical expression. 

For studying the variation of the storage delay time and fall time with 
respect to lifetime (jn and Tp), we have taken a long base pn diode with 
cm"®, so that we can see the minority carrier distribution prop- 
erly (in case of n'^p diode, holes on the n side are very small as compared to 
electrons on the p side) on the falling edge of the input voltage. Figure 3.30 
shows the decay of minority carriers for different values of r„ and Tp. In Fig 
3.30a, Tn = Tp, we see that on both sides decay of excess minority carriers is 
occurring at the same rate and reaching zero at the same time. In Fig 3.30b, 
Tn = 4 Tp, hence electrons on p side are decaying slowly compared to holes 
on n side. In Fig 3.30c, rp= 4 t„, therefore the holes on n side decay slowly 
compared to electrons on p side. 

The peak concentration of the minority carriers at the edge of the depletion 
region boundary is also increased. It can be explained as follows: suppose the 
forward current is assumed to be constant at ~ (VsIRl), then Jn = 
and from this equation we can write 6n{x) oc Therefore, the excess 

minority carriers will increase as r„ is increased, while the current is constant. 
A similar explanation for Fig 3.30c, where Tp = 4 Tn, can be given. The 
variation of the storage delay time and fall time with respect to r„/rp is shown 
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in Figs 3.31 and 3.32 respectively. 


3.3 Comparison with SPICE 

For comparison with SPICE [11], a short base diode is taken for simulation 
because SPICE uses short base diode model. The circuit used is shown in 
Fig 2.25c ( where forward voltage= 10 ^blts, reverse voltage= -2 Volts). The 
parameters required in SPICE, i.e., transit time (TT = where Wj, is the 
width of p side, junction voltage (VJ), reverse saturation current (IS), parasitic 
resistance (RS), and zero bias pn capacitance (CJO) are kept the same as those 
which are used in simulation. The current versus time plot is shown in Fig 

3.33. 

It is seen that SPICE gives very high value of the storage time delay. It 
is almost five times the simulated value. On the other hand SPICE gives 
negligible fall time while simulated results gives finite amount of fall time . 
This difference is probably due to the fact that, the decay of minority carrier 
concentration as a function of time is rather complicated as shown if the Fig 

3.34. In SPICE, only simple approximation is used to describe the minority 
carrier concentration. The difference for fall time delay is due to the improper 
capacitance modelling. 
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Figure 3.22: (a) Input Volta^ vs time; (b) Diode Voltage vs time; (c) Diode 
Current vs time; (d) Electric field nem the junction vs time. 






Figure 3.23: Turn off transient of the diode current vs time for different exter- 
nal load, keeping the forward and the reverse bias constant. 



Figure 3.24: Turn off transient of the diode current vs time tor different reverse 
biases, keeping the forward bias constant. 
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Figure 3.25: Storage time delay vs Reverse bias. 



Figure 3.26: Fall time delay vs Reverse bias. 
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Figure 3.29: Fall time delay vs Forward bias. 
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(c) 

Figure 3.30: Minority carriers concentration vs distance for different values of 
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time delay vs Tnjrp. 



Figure 3 . 32 ; Fall time delay vs r„/rp. 


39 





Figure 3.33: Diode current vs time. 



after the negative edge of the voltage step. 
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Chapter 4 


Conclusion and scope for future 
work 


In the present work, numerical solution of basic semiconductor equations un- 
der transient conditions is presented. Carrier concentrations and potential in 
the interior of the device for time dependent current or voltage excitation are 
shown. This analysis allows a thorough understanding of the basic time de- 
pendent phenomena governing the physical behaviour of the semiconductor 
device under typical external excitations. 

We have found good agreement between our simulation results and theory 
with respect to storage delay time. On the other hand, analytical expression 
was found to underestimate the fall time delay substantially. SPICE sunulation 
gives very large storage delay time while it does not gives any fall time delay. 
The steady state solution and the transient solution represent a powerful tool 
of investigation in areas of computer aided device modelling, in 'particular the 
analysis of devices under switching conditions. 



9 

In the present work, complications such as Auger recombination, nonuni- 
form doping density, etc., are not considered because we wanted to compare our 
results with the analytical expressions w'hich are also derived under simple as- 
sumptions. To come close to the real diode more extensive studies are required 
so that all the effects are considered. As the numerical simulation results are 
more accurate, improved models for practical diodes can be suggested. 
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